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Abstract 

The dominant contributions to W^H"^ production at the LHC are the tree-level 
bb annihilation and the gg fusion. We perform for the case of the complex MSSM 
a complete calculation of the NLO EW corrections to the bb annihilation channel 
and a consistent combination with other contributions including the standard and 
SUSY QCD corrections and the gg fusion, with resummation of the leading radiative 
corrections to the bottom-Higgs couplings and the neutral Higgs-boson propagators. 
We observe a large CP-violating asymmetry, arising mainly from the gg channel. 



1 Introduction 



The discovery of charged Higgs bosons at the running Large Hadron Colhder (LHC) would 
be unambiguous evidence for new physics. Important mechanisms to produce include 
gh — )■ tH~] qq,gg — ?■ H^H~ and hh^gg — )■ W^H^ (see [l] for a review and references). 
The first and the last channel are of particular interest since they allow to search for CP- 
violating effects at the LHC associated with physics beyond the Standard Model. Recently, 
the CP- violating asymmetry for tH~ /iH^ production has been calculated in [2]. 

There have been many discussions devoted to the pp — )■ W^H^ processes in the 
Minimal Supersymmetric Standard Model (MSSM) over the last two decades. These 
studies assume that all the soft supersymmetry-breaking parameters are real and hence 
CP violation is absent. The two main partonic processes are bb annihilation and the loop- 
induced gg fusion. The first study js] computed the tree-level bb contribution and the 
gg process with third-generation quarks in the loops using rrif, = approximation. This 
calculation was then extended for finite mi,, thus allowing the investigation of the process 
for arbitrary values of tan (3 (the ratio V2/V1 of the vacuum expectation values of the two 
Higgs doublets) (4||5]. The inclusion of the squark-loop contribution to the gg channel was 
done in (6||7]. The next-to-leading order (NLO) corrections to the bb annihilation are more 
complicated and not complete as yet; the full NLO electroweak (EW) corrections are still 
missing. The Standard Model QCD (SM-QCD) corrections were calculated in l8l|9), the 



supersymmetric-QCD (SUSY-QCD) corrections in 10,11 , and the Yukawa part of the 



electroweak corrections in 12 . There are also studies on the experimental possibility of 



observing VV^H^ production at the LHC with subsequent hadronic H tb decay 13 



and leptonic H — )• r z/^ decay 14 , 15 



The aim of this paper is multifold. First, we extend the calculation for pp — )■ W^H^ to 
the MSSM with complex parameters (complex MSSM, or cMSSM). Second, the full NLO 
EW corrections to the bb annihilation channel are calculated and consistently combined 
with the other contributions to provide the complete NLO corrections to the pp — )■ W^H^ 
processes. Third, CP-violating effects arising in the cMSSM are discussed. The important 
issues related to the neutral Higgs mixing and large radiative corrections to the bottom- 
Higgs couplings are also systematically addressed. 

In the cMSSM, new sources of CP violation are associated with the phases of soft- 
breaking parameters and of the Higgsino-mass parameter /i. Through loop contributions, 
CP violation also enters the Higgs sector, which is CP conserving at lowest order (see for 
example 16 for more details and references). As a consequence, the /i, H and A neutral 
Higgs bosons in general mix and form the mass eigenstates /ii,2,3 with both CP even and 
odd properties, which can have important impact on many physical observables. 

The bottom-Higgs Yukawa couplings are subject to large quantum corrections in the 
MSSM. We use the usual QCD running bottom-quark mass to absorb large QCD correc- 
tions to the LO results. The potentially large SUSY-QCD corrections, in the large tan (3 
limit, are included into the quantity Amj,, which is complex in the cMSSM and can be 
resummed (Section 2.1). 

The paper is organized as follows. Section[2]is devoted to the subprocess bb — )■ W^H^ , 
including the issues of effective bottom-Higgs couplings and neutral Higgs mixing. The 
calculation of the gg fusion part is shown in Section [3j Hadronic cross sections and CP- 



1 



violating asymmetry are defined in Section |4} Numerical results are presented in Section [5] 
and conclusions in Section |6| Feynman diagrams, counterterms, and renormalization 
constants can be found in the Appendices. 



2 The subprocess bb W^H 



± 




Figure 1: Tree-level diagrams for the partonic process bb — )■ W^H^. hi with i 
denote the neutral Higgs bosons h, H and A, respectively. 



1,2,3 



At the tree level, there are four Feynman diagrams including three s-channel diagrams 
with a neutral Higgs exchange and a t-channel diagram, as shown in Fig. [TJ The tree-level 
bottom-Higgs couplings read as follows. 
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where Pl,r = (1 =F75)/2, sw = sm9w, and a is the tree-level mixing angle of the two CP- 
even Higgs bosons. In order to obtain reliable predictions, two important issues related to 
the bottom-Higgs Yukawa couplings and the neutral Higgs mixing have to be addressed. 
These quantities can get large radiative corrections as will be detailed in the next two 
sections. 



2.1 Bottom-Higgs couplings 

In the context of the MSSM, the bottom-Higgs couplings can get large SM-QCD, SUSY- 
QCD and EW corrections. These large universal corrections can be absorbed into the 
bottom-Higgs couplings in two steps. First, the SM-QCD corrections are absorbed by 
using the running bottom-quark mass at one-loop order via 
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We note, in passing, that the relation between the pole mass and the MS mass is different 
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It can be proved that by using the running bottom-quark mass in Eq. ([T]) the SM-QCD one- 
loop corrections are independent of In(m^) 17 . We will therefore replace = mP^(/i) 
Eq. ^ 
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can be related to the QCD-MS mass mb{mh), which is extracted from 
experimental data and is usually taken as an input parameter, at two-loop order as follows 
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where n is the number of active quark flavours and the MS running mass is evaluated 
with the two-loop formula 



U(i{fiB.,mt)U5{mt,mb)mb{mb) for fin > nit 
U5{fiR,mb)fnb{fnb) for fiR < rrit 

where the evolution factor [/„ reads (see e.g. [I9] ) 
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The second step is to absorb large universal SUSY-QCD and EW corrections into 
the couplings in Eq. Q. This is achieved by using the following effective bottom-Higgs 
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Al 



1 — Ab/(tan /3tan a) 

1 + A, 
1 + Afotan a / tan /3 

l + Ab ' 
1 - Af,/(tan/3)2 

1 + Ab ■ 
3 



The leading corrections proportional to C(astan (3, attan (3, atan (3) with at = hf/{4:7i) and 
ht being the superpotential top coupling are included in Arub (l9). This quantity is UV 
finite and can be calculated by considering the one-loop corrections to the Hi^bb coupling 
(which is zero at tree level) where H2 is the neutral component of the second Higgs 



doublet. It can also be extracted from the one-loop bottom-quark self-energy 23,24 . In 
the cMSSM, we find 

Am, = Amf^^ + Amf^, 
^^SQCD ^ ^M^MgV tan/3 /(m?^,m?^,m|), Q = {m-,^ + m-,^ + m^)/3, 

Amf^^ = Amf * + Amf + Amf , 

Amf = ^A>nan/3 J(m|,m|,|/xn 

^^f = -^M*/i*tan/3[2|t/fi|2j(m|, \M,\\ + 2|f/*j2/(^|, \m,\\ \f,\') 

+ \ul\'l{ml, l/ip) + \ul\'l{ml, |/i| 



+3 {2\U',,\' + |f/2\r)/K, |Mi|2, + 2/(m[, m^^, \M^\')] , (9) 
with the auxiliary function 

1 f a b c\ 

I(a, b, c) = —7 —7- —. r I a6 In - + 6c In - + ca In - . (10) 

[a — b)[b — c)[c — a) \ b c a J 

Ml, M2, M3 (each with a phase Mj = \Mj\e'^'^^) and /i = |/i|e*''^'" are the bino (B), wino 
(W), gluino (g) and Higgsino (H) mass parameters, respectively. Af = \Af\e^'^f, here / 
means fermion, denotes the soft supersymmetry-breaking trilinear scalar coupling, bi and 
ii with i = 1,2 are the sbottom and stop mass eigenstates, respectively. U'^ and are 
2x2 mixing matrices. By setting all the phases to zero we obtain the results for the real 



MSSM (rMSSM), which agree with those given in 19,25 . Since we are also interested in 



the effect of the Ab phase, corrections proportional to Ab are resummed by 20 , 26 

, Am, 

= ITaI' 

Ai = -^-^M;AbIiml,mlml). (11) 

We remark that A, is complex and depends on (^^, (pf, (pi with i = 1,2,3. The effective 
couplings in Eq. ([T]) are used in the calculations of the tree-level, SM-QCD and SUSY- 
QCD contributions to the bb — )■ W^H^ process and the gg fusion. For the NLO EW 
corrections we use the tree-level couplings Eq. ([l]) with nib = 'mf^il^R)- 

In the explicit one-loop calculations, we have to subtract the Afe-related corrections 
which have already included into the tree-level contribution to avoid double counting. 
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This can be done by adding the following counterterms 
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to 5mfe in the corresponding bottom-Higgs-coupling counterterms, as listed in Appendix [B| 
Moreover, Eq. ^ is used with A^ = Amf^^^, Amf^^ for the SUSY-QCD and EW 
corrections, respectively. 

2.2 Neutral Higgs-boson propagators 

In the MSSM, the neutral Higgs boson masses are subject to large radiative corrections 
in particular from the Yukawa sector of the theory. As a consequence, the tree-level Higgs 
masses can be quite different from the physical ones. This important effect should be 
considered in the NLO calculations of processes with intermediate neutral Higgs exchange. 

In our calculation, both subprocesses include s-channel diagrams with internal neutral 
Higgs bosons, (Fig. [l] and Fig. 20). For hh — > W^H^ there is also a t-channel diagram 
at tree level which gives the dominant contribution at high energies. Thus, the higher- 
order corrections to the internal Higgs propagators are not expected to have important 
effects in this subprocess at high energies. This will be verified in our numerical studies 



in Section 5.3 The situation is different with the gg fusion since the s-channel (triangle) 



contribution is large. The higher-order corrections to the internal Higgs propagators can 



therefore be significant in this will be confirmed in Section 15.51 This issue has 

not been addressed in the previous studies. 

In a general amplitude with internal neutral Higgs bosons that do not appear inside 
loops, the structure describing the Higgs-exchange part of an amplitude is given by 



^(/) = ^r,A,,(/)r,-, i = h,H,A, 



(13) 



where Fjj are one-particle irreducible Higgs vertices, p is the momentum in the Higgs 
propagator, which is given in terms of the 3x3 propagator matrix 



A(p2) = ^[p2_M(p2)]-i, 

'ml - S,,(p2) -S,^(p2) -tUP^) 

M(/) = I -t^H{p') ml-tHH{p') -tuAip') 

-thAip") -tHAip") m\ - tAAip") 
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rrii 



h, H, A) are the lowest-order Higgs-boson masses, and Sij the renormalized 



self-energies. The physical masses can be found by diagonalizing the above matrix 27 



By using this propagator matrix we effectively resum all the one-loop corrections to the 
neutral Higgs self-energies. 

In our calculation, we keep the full propagator matrix and Eq. ( 13 ) whenever neutral 
Higgs-bosons are exchanged connecting three-point vertices in the tree-level bb contri- 
butions and in the gg fusion diagrams. As a consequence, when including the NLO EW 
corrections, we have to discard all Feynman diagrams containing diagonal and nondiagonal 
hihj self-energies to avoid double counting (see Fig. 19). Whenever neutral Higgs bosons 
appear inside a loop, the tree-level expressions are used for propagators and couplings. 



The renormalized Higgs self-energies in Eq. (14) are calculated at NLO by using the 
hybrid on-shell and DR scheme (see Section 2.6 and 27 for details). Our results have been 



successfully checked against the ones of FeynHiggs |27f 30] . It is noted that FeynHiggs has 
the option to include the leading two-loop 0{asat) corrections in the cMSSM 31,321. We 



have verified that the effects of these two-loop corrections are negligible in our numerical 
analysis and we thus chose to perform the numerical evaluation with the one-loop self- 
energies. 

To quantify the effect of the neutral Higgs propagators we introduce two approxima- 
tions for the subprocess bb — )■ W^H^: The improved-Born approximation (IBA) including 
both the Afe resummation and the neutral Higgs mixing resummation, and the simpler 
version IBAl which contains only the resummed A;, together with tree-level Higgs boson 
masses and couplings. By LO we refer to the tree-level bb W^H^ contribution with 



m 



DR 



(Hr) and the tree-level Higgs sector. 



2.3 SM-QCD corrections 



The NLO contribution includes the virtual and real gluonic corrections. The virtual cor- 



rections, displayed by the Feynman graphs in Fig. 16, contain an extra gluon in the loops. 
The calculation is done by using the technique of constrained differential renormalisation 
(CDR) 33 which is, at one-loop level, equivalent to regularization by dimensional reduc- 
We have also checked by explicit calculations that it is also equivalent to 



35 



tion |34 

dimensional regularization ||36j in this case. 

Concerning renormalization, the bottom-quark mass appearing in the Yukawa cou- 
plings is renormalized by using the DR scheme. It means that the running m^^lfipd (see 



Section 2.1) is used in the Yukawa couplings and the one-loop counterterm reads 



6m 



DR 



Ait 



(15) 



where Cp = 4/3, Cuv = ~ 7^ + ln(47r) in D = 4 — 2e space-time dimensions with 
'jE denoting Euler's constant. The bottom-quark mass related to the initial state (in 



the kinematics 



mt 



and the spinors) is treated as the pole mass since the correct 



on-shell (OS) behavior must be assured. Indeed the m^^ effect here is very small and can 
be neglected. As mentioned in Section 2.1, the final results are independent of In(mP^). 



We will therefore set = m^^lfiR) everywhere in this paper. The finite wave- function 
normalization factors for the bottom quarks can be taken care of by using the OS scheme 
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for the wave- function renormalization. For the top quark, the pole mass is used throughout 
this paper. Accordingly, the mass counterterm is calculated by using the OS scheme 
(Appendix [B]). 

The real QCD corrections consist of the processes with external gluons, 

b + b W-+H+ + g, 
b + g b + H+ + W-, 

b + g b + W- + H+, (16) 



corresponding to the Feynman diagrams shown in Fig. 17 For the gluon-radiation process, 
soft and coUinear divergences occur. The soft singularities cancel against those from 
the virtual corrections, while the coUinear singularities are regularized by the bottom- 
quark mass. The gluon-bottom-induced processes are infrared finite but contain coUinear 
singularities, which are regularized by the bottom-quark mass as well. After adding the 
virtual and real corrections, the result is coUinear divergent and proportional to ln(m^/s), 
where a/S is the center-of-mass energy. These singularities are absorbed into the bottom 
and gluon parton distribution functions (PDF), as discussed in Section 111 



Following the line of 37 , we apply both the dipole subtraction scheme 38 , 39 and 



the two-cutoff phase space slicing method 40 to extract the singularities from the real 



corrections. The two techniques give the same results within the integration errors. How- 
ever, the error of the dipole subtraction scheme is much smaller than the one of the phase 
space slicing method. We will therefore use the dipole subtraction scheme in the numerical 
analysis. 



2.4 Subtracting the on-shell top-quark contribution 



A special feature of the gluon-induced processes in ( 16 ) is the appearance of on-shell top- 



quarks decaying into bW (and bH~^ when kinematically allowed), which requires a careful 



treatment and has been discussed in the previous literature, e.g. in 41 43 . Our approach 



is similar to the one described in 41 , 42 , with the difference that we perform the zero 
top-quark width limit. 

We demonstrate the procedure in terms of the process bg — > W^H^b. The Feynman 
diagrams (Fig. 17:) include a subclass involving the decay i — )■ bW~ . When the internal 
t can be on-shell, the propagator pole must contain a finite width F^, which is regarded 
here as a regulator: 



— rrif + irritTt 



(17) 



This on-shell contribution is primarily a iH^ production and should therefore not be 
considered a NLO contribution. For the genuine NLO correction, the on-shell top contri- 
bution has to be discarded in a gauge invariant way. Starting from the full set of diagrams, 
the squared matrix element reads as follows. 



|Mp = iMosr + 2Re[MosM,Vos] + |M,on-os| 
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where the subscripts os and non-os denote the contribution of the on-shell i diagrams and 
the remainder, respectively. The OS part, differential in the bW invariant mass, to be 
subtracted can be identified as 



• W-H+b 



sub p 

= >-^H+^Br(t->&Py-) (19) 



where Br(t — )■ hW~) = F^^g^./Ft. The ratio on the right-hand side (rhs) of Eq. (19) 
approaches 6{M^^^ — ml) when Ff — )■ 0. The subtracted NLO contribution, regularised 
with the help of F^, can be written in the following way, 

J>g^w-H+-b(j.~. _ f ...2 ( daps^ -b,^H+t rntTtBTjt bW-) 

+ -^"""''+-^:gr"'^ (20) 



where the interference and non-OS terms arise from the second and third terms in Eq. ( 18 ). 



There is strong cancellation between the first term in the rhs of Eq. (20) and the rest after 
subtraction of the collinear part, which makes the result of Eq. (20) very small, yielding 
an essentially linear dependence on F^ as displayed in Fig. |2} We can thus perform the 
limit Fj — 7- and obtain a gauge invariant expression by 

^h->w-Hn ^ aJ,r^-^"^(F,). (21) 
^ rt-5>o ^ 

Fig. |3] shows that the finite gluon-induced contribution obtained in this way at the 
hadronic level (after proper subtraction of the collinear part) is very small for large values 
of Mh±, but it can be of some significance when the charged Higgs boson is light. 

The method described above is completely analogous for the process bg — )■ W~H^b. 
For low masses, Mh± < mt, the intermediate on-shell top quark can also decay into 
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Figure 3: The finite hadronic cross section cr\l'^^ after subtracting the OS top-quark 
and the collinear- singularity contributions as a function of Mh± ■ 



H^h. This additional OS contribution can be extracted by using the same extrapolation 
method. For completeness, we list here the expressions for the decay widths of t — > hW^ 
and t — 7- hH^ at lowest order, 
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where the 6-quark mass has been neglected 
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2.5 SUSY-QCD corrections 

The NLO SUSY-QCD contribution consists only of the virtual one-loop corrections, visu- 



alized by the Feynman diagrams with gluino loops in Fig. [T8j The only divergent part is 
the top-quark self energy, which is renormalized in the on-shell scheme. As discussed in 



Section 2.1, large corrections proportional to asM|/i*tan/3 have been summed up to all 
orders in the bottom-Higgs couplings included in the IBA. We therefore have to subtract 
this part from the explicit one-loop SUSY-QCD corrections to avoid double counting. 



2.6 Electroweak corrections 

The full NLO EW contributions to the processes hb — > W^H^ in the cMSSM have not 
been computed yet. They comprise both virtual and real corrections. For the virtual 



part. Fig. 19 illustrates the various classes of one-loop Feynman diagrams. As before, 
the calculation is performed using the CDR technique. We have also worked out all the 



necessary counterterms in the cMSSM and implemented them in FeynArts-3.4 44,45 



Explicit expressions for the counterterms can be found in Appendix |Bj For the Higgs field 
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renormalization and tan /3, we use the DR renormalization scheme as specified in 27 . 



Hence, the correct OS behavior of the external must be ensured by including the 



finite wave-function renormalization factor 46 



(24) 



where Tjh-h+ {p^) is the renormalized self-energy. The other renormalization constants 
are determined according to the OS scheme. To make the EW corrections independent of 
Inmj from the fight fermions / 7^ t, we use the fine-structure constant at M^, a = a{Mz) 
as an input parameter. This means that we have to modify the counterterm as 



Aa(M| 
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~ 2 

fc2=0 



A«(M|), 
_ ReS^^(M|) 



(25) 



where the photon self-energy includes only the light fermion contribution, to avoid double 
counting. 

The real EW contributions correspond to the processes with external photons. 



b + b 
b + -f 
b + -f 



— > 



W- +H+ + 7, 
b + H+ + W-, 
b + W- + 



(26) 



described by the Feynman diagrams of Fig. 21 They are calculated in the same way as the 



real QCD corrections, discussed in Section 2.3 and Section 2.4 Naively, we would expect 
this photon contribution to be much smaller than the one from the gluon, due to the 
smallness of the EW coupling a and the photon PDF. This is not always true, however, 
since the photon couples to the and as well. The soft singularities are completely 
cancelled, as in the case of QCD. The EW splitting 7 — )■ H^H~ (similarly for 7 — t- 
W^W~), on the other side, can introduce large collinear correction in the limit Mu± /Q — ?■ 
0, Q is a typical energy scale. The constraint > Mw prevents those splittings from 
becoming divergent. We observe, however, that the finite corrections (after subtracting 
the collinear bottom-photon and the OS top-quark contributions) from the above 67 
process are still larger than the corresponding QCD ones for Mh± < 200 GeV, e.g. for 
Mh = 150 GeV and a/s = 14 TeV by a factor of 2. The photon-induced contribution 
should thus be included in the NLO calculations for /H^ production at high energies. 
This requires the knowledge of the photon density in the proton, which at present is 



contained in the set MRST2004qed 47 of PDFs. 



3 The subprocess gg W^H^ 

The subprocess gg — )■ W^H^ is loop induced, in the MSSM with quark- and squark-loop 



contributions. Fig. 20 summarizes the various one-loop Feynman diagrams, which involve 
three- and four-point vertex functions. Since the (s) quarks are always coupled to a Higgs 
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boson, the one- loop amplitude is proportional to (s)quark-Higgs couplings. The dominant 
contributions therefore arise from the diagrams with the third-generation (s) quarks. As 
in [7], the contribution from the first two generations of (s)quarks is neglected in this paper. 
Compared to the previous work our calculation is improved by using the effective 
bottom-Higgs couplings and the resummed neutral Higgs propagators. It turns out that 
these improvements sizably affect both the cross section and CP-violating asymmetry. 
We have checked our results against those of [7j for the case of the real MSSM using the 
tree- level couplings and Higgs propagators and found good agreement. 




Figure 4: Feynman diagrams that can produce three-point Landau singularities. 



We notice an interesting feature related to the anomalous thresholds. Fig. lb of fl\ 
shows a very sharp peak close to the normal ti threshold. Careful observation reveals that 
the peak position is slightly above 2mt and is obviously more singular than the normal 
thresholds in Fig. la of |7|. This is indeed an anomalous threshold corresponding to the 



three-point Landau singularity (see 48 ,49 and references therein) of the triangle and box 



diagrams in Fig. |4j A simple calculation following [48] yields the peak position at 

Speak = ^[{Ml^+Mlr){ml + ml)-{ml-m'lf -Ml^Ml, 

-X^ml ml MlJX^ml ml M^)l , (27) 



with A(x, y, z) = + y'^ + — 2{xy + yz + xz). The partonic cross section is divergent at 
s = Speak but the result is finite at the hadronic level, i.e. after integrating over s, since 
this singularity is logarithmic and thus integrable. The conditions for this anomalous 



threshold to be in the physical region can also be given 48 



^' nib 



mb + nit < Mh± < \j2{nii + ml) - Ml,. (28) 
Similarly, the three-point Landau singularities can occur in the squark diagrams. 

4 Hadronic cross section and CP asymmetry 

The LO hadronic cross section, in terms of the LO partonic hh annihilation cross section, 
is given by 

= j dx^dx2[F^{x,,fiF)F-^{x2,fiF)afoic^lfiR) + (1 ^ 2)], (29) 
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where F^^-^{x, hf) is the bottom PDF at momentum faction x and factorization scale /i^. 
Other gg-subprocesses (g = u, d, c, s) are neglected due to the smallness of light-quark- 
Higgs couplings. 

The NLO hadronic cross section reads as follows, 



a 



pp 

NLO 



E — 



+ (1^2)] 



(30) 



where i,j = {b,b,g,-y) and 



^NLO = ^IBAi.(^ ) + AsM-QCD<3-^Lo(« ^^s) + AsuSY-QCDO-^lo (« 



(31) 



contain the various NLO contributions at the parton level, discussed in the previous 
sections. As mentioned there, the mass singularities of the type as ln(mf,) and a ln(m6) 
are absorbed in the quark distributions. We use the MRST2004qed set of PDFs [47J, which 
include 0{as) QCD and 0{a) photonic corrections. As explained in 50 , the consistent 
use of these PDFs requires the MS factorization scheme for the QCD, but the DIS scheme 
for the photonic corrections. We therefore redefine the (anti-)bottom PDF as follows, 



q{x) 



qix,f4) 



O^sCf 
271 



^ dz (X 



z \z 



2 

, l^F 



In 



[P,,(^)(ln(l-z)2 + l)]+ + C,7(.) 



olQ\ dz (X 
— g - 
z \z 



27r 



2 

, l^F 



In 



[P,,(^)(ln(l-z)2 + l)]+ + C, 



DIS/ 
qq \ 



27r 
3aQ| 
27r 



^ dz (X 

—9\- 

z \z 

^ dz (X 

— 7 - 
z \z 



2 

) l^F 



2 

, l^F 



In 



In 



Pqg + ^^^(^) 



■yDIS/ 



(32) 



with Cf = 4/3, Tp = 1/2. The splitting functions are given by 

1 + ^2 



Pqq{^) 



1 - z 



(33) 



and the [...]+ prescription is understood in the usual way, 



dzf{z) 



9{z) 
1 - z 



',J/M-/(i)bM_^(,) 

1 ^ 



dz 



1 - z' 



(34) 
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Following the standard conventions of QCD, the factorization schemes are specified by 



ln(- 



^pis, 



<?7 



Z) 



I- 



V-z 3 
~ 4 

- 8^^ + 8z 



9 + 5^ 



(35) 



Having constructed in this way the hadronic cross sections a{pp — )■ W^H^), we can define 
the CP-violating asymmetry at the hadronic level in the following way, 



^CP 

pp 



a{pp W-H+) + a{pp W+H- 



(36) 



The numerator gets contributions from the NLO-66 corrections (the LO is CP conserving) 
and the loop-induced gg process. However, the latter is much larger than the former due 
to the dominant gluon PDF. The CP- violating effect is therefore mainly generated by the 
gg channel. The LO-bb contribution adds only to the CP invariant part and therefore 
reduces the magnitude of the CP asymmetry. 



5 Numerical studies 



5.1 Input parameters 



We use the following set of input parameters for the SM sector 51 , 52 
asiMz 



Mw 
mt 



0.1197, 
80.398 GeV, 
173.1 GeV, 



a{Mz) 
Mz 
mb(mb) 



1/128.926, 
91.1876 GeV, 
4.2 GeV. 



(37) 



a 



MS 



{(ir) at three-loop order 



51 . mbirrib) is the QCD-MS 6-quark 



We take here as 

mass, while the top-quark mass is understood as the pole mass. CKM matrix elements 
are approximated by Vtd = Vt^ = and Vtb = 1. 

For the soft SUSY-breaking parameters, we use the adapted CP-violating benchmark 



scenario (CPX) p,54 



1^1 = 2 TeV, IM2I = 200 GeV, IM3I = 1 TeV, \At 



I A 



\Ar\= 900 GeV, 



M^ = M^ = Mr, = M? = Me. = Mc 



SUSY 



500 GeV. 



(3^ 



Since the Yukawa couplings of the first two fermion generations proportional to the small 
fermion masses are neglected in our calculations, we set A/ = for / = e, /i, m, (i, c, s. The 
values of Mi and M2 are connected via the GUT relation \Mi\ = 5/3 tan^ IM2I. We 
can set 02 = while keeping 0i as a free parameter. The complex phases of the trilinear 
couplings At, Ah, A^- and the gaugino-mass parameters Mj with i = 1,2,3 are chosen as 
default according to 



TT 

2' 



(39) 
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unless specified otherwise. The phase of /i is chosen to be zero in order to be consistent 
with the experimental data of the electric dipole moment. We will study the dependence 
of our results on tan /3, Mh±, (pt and 03 in the numerical analysis. The (pt, dependence is 
not very interesting since it is similar to but much weaker than that of (pf 

The scale of in the SUSY-QCD resummation of the effective bottom-Higgs couplings 
Eq. (joj) is set to be Q = (m^^ + m^^ + mg)/3. If not otherwise specified, we set the 
renormalization scale equal to the factorization scale, fiR = fip, in all numerical results. 
Our default choice for the factorization scale is fipo = My/ + Mh±. 

Our study is done for the LHC at 7 TeV and 14 TeV center-of-mass energy. In the 
numerical analysis, we will focus on the latter since the total cross section is about an 
order of magnitude larger. Important results will be shown for both energies. 



5.2 Checks on the results 



The results in this paper have been obtained by two independent calculations. We have 

fferent Fortran 77 
The phase-space 



produced, with the help of FeynArts-3 . 4 and ForniCalc-6 . 35), two different Fortran 77 
codes. 



Loop integrals are calculated by using LoopTools/FF 35,55 



integration is done by using the Monte Carlo integrators BASES 56 and VEGAS 57 . The 



results of the two codes are in full agreement. On top, we have also performed a number 
of other checks: 

For the process gg — )■ W^H^, we have verified that the results are QCD gauge in- 
variant. This can be easily done in practice by changing the numerical value of the gluon 
polarization vector e^(]9, g), where p is the gluon momentum and q is an arbitrary refer- 
ence vector. QCD gauge invariance means that the squared amphtudes are independent 
of q. More details can be found in 
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As already mentioned, we compared our results 
also to the ones of j^for the rMSSM and obtained good agreement. 

For the process hh — )■ W^H^, besides the common checks of UV and IR finiteness, 
we compared our virtual EW corrections to those obtained by using SloopS 59 ,60 , and 



the SUSY-QCD corrections to the results of Ranch 11 for the case of vanishing phases. 
Again, good agreement was found. 



5.3 pp /hh W^H^: LO and improved-Born approximations 

In this section, we study the effect of the bottom-Higgs coupling resummation described 



in Section 2.1 and of the Higgs propagator matrix discussed in Section 2.2 

The results for the approximations IBA and IBAl defined in section 2^ are illustrated 
in Fig. [5] showing the dependence on tan (5 in the left panel and on the mass Mh± in the 
right panel. The relative correction 5, with respect to the LO cross section, is defined as 
5 = (c"iBA — crLo)/<^LO- Fo^ Small values of tan /3 the left-chirality contribution proportional 
to rrit/ tan/3 is dominant while the right-chirality contribution proportional to mf,tan/3 
dominates at large tan /3. The cross section has a minimum around tan /3 = 8. 

The effect of Ab resummation is best understood in terms of Fig. [5] and Fig. [6j The 
important point is that A;, is a complex number and only its real part can interfere 
with the LO amplitude. Thus, the A;, effect is minimum at cpt.z = ±7r/2 where the 
dominant Am^*^'"'^'^* are purely imaginary and is largest at (j)t,3 = 0, ivr. 0j enters via 
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Figure 5: The leading order (LO) cross section with nif, = m^^- and the two improved 
Born approximations (IBA) as functions o/tan/3 (left) and Mh± (right), ctibai includes 
the Ah resummation but not the Higgs mixing resummation, while (Jiba includes both. 
The lower panels show the corresponding relative corrections with respect to the LO result. 
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Figure 6: Similar to Fig. but with (pt (lefi) and 03 (right) varied instead. 



15 




'2 4 6 8 10 12 14 16 18 20 100 150 200 250 300 350 400 450 500 550 600 

tanp M^^, [GeV] 



Figure 7: The cross section obtained by using IBA and including various nonuniversal 
NLO corrections as functions of tan (3 (left) and Mh± (right). The lower panels show the 
corresponding relative corrections to the IBA result. 



EW corrections and 03 via the SUSY-QCD contributions. Fig. [6] shows that the A;, effect 
can be more than 150%. In Fig. [s] where is mostly imaginary we see the effect of 
order 0{Al) which is about —15% at tan/3 = 10. We also observe that the Higgs mixing 
resummation in the s-channel diagrams has a much smaller impact, less than 10%, as 
expected. 

5.4 pp/bb VF^ff^: full NLO results 

In this section, we investigate the effects of the SUSY-QCD, SM-QCD, and EW contri- 
butions at NLO. As in the previous section, we present here two sets of plots. In Fig. [7] 
we show the dependence of the total cross sections on tan /3 and Mh± at the default CPX 
phases, in particular 4>t = 4>3 = 7r/2. As explained above, the 0{Ab) effect is turned 
off in this CPX scenario. The SUSY-QCD and EW NLO terms are therefore small at 
large tan /3, as shown in Fig. [T] (left). The SM-QCD correction is about —20% for small 
tan P and changes the sign around tan /3 = 11 due to the competition between the bb and 
the gf-induced contributions. All the NLO contributions for different values of tan /3 and 
Mh± can be found in Table [T] Fig. [8] shows the dependence of the total cross sections 
on (pt and 03 for tan /3 = 10 and Mh± = 200GeV. The EW corrections depend strongly 
on 04, and the SUSY-QCD corrections on 03. At 4>t = (ps = 0, ±7r the effects are largest. 
The remaining EW and SUSY-QCD corrections, beyond the 0{Ai,) contribution, are still 
rather large. In particular, there is the following term of the SUSY-QCD correction. 
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Figure 8: Similar to Fig. ^ but with (pt (left) and 03 (right) varied instead. 




Figure 9: Diagrams that can introduce large SUSY-QCD (left) and EW (right) corrections, 
are the charged Goldstone bosons. 
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g^-M3>nan/3J(m^), 
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which can be included in the top- Yukawa part of charged Higgs couphngs as follows 

ie ( nit 



(40) 



KbH- 



\p2swMw vtan/3 
ie 



;i-A0PL + <«tan/3AfP, 



w 



m°^tan/3A3p^+ 



tan (3 



(41) 



This term originates from the left diagram in Fig. [9] and is important for small tan/3. 
This finding agrees with the discussion in ^26j where other subleading corrections are also 
discussed. If the couplings Eq. (41) are used we find that the new-improved LO results 



move significantly closer to the full NLO results in Fig. |8] (right). The situation in the 
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Figure 10: The cross section (left) and CP asymmetry (right) as functions of i 



>t- 



left part of Fig. [8] is due to the EW corrections. It indicates that there are still large 
corrections proportional to Atfiat/ (4:71) which can be associated with the right diagram 
in Fig. [9j 

The SM-QCD corrections (and EW corrections to a lesser extend) have a striking 
structure for small masses Mh± < rrit (Fig. [7| right part). This is due to the finite 
contribution of the process hg — )■ H^h. When Mh± < rrit the intermediate top quark 
can be on-shell and can decay to H^b. As discussed in Section 2^, this OS contribution 
has to be properly subtracted. The structure indicates that the OS top-quark effect 
cannot be completely removed and this quantum effect on the W~H~^ production rate is 
an interesting feature, which was not discussed in previous studies [SlIO]. 



5.5 pp/gg W^H^: neutral Higgs-propagator effects 

Even though the gg-fusion subprocess is loop induced, its contribution can be of the same 
order as the tree-level bb — )■ W^H^ contribution. Neutral Higgs bosons are exchanged 
in the s-channel and can be described by using effective bottom-Higgs couplings and the 
full Higgs-propagator matrix. The impact of the latter on the total cross section and CP 
asymmetry is large as can be seen from Fig. 

20% at (pt = ±7r, while the CP asymmetry increases about 25% at 



The cross section can be reduced by 
= ±7r/2. This is 

consistent with the discussion in Section 2^ We also observe that the gg contribution is 
very sensitive to (pf 
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Table 1: The total cross section in fb for pp/bb -^■ W^H^ including the IBA and various 
nonuniversal NLO corrections and for pp/ gg — )• W~H^ at y/s = 14 TeV. The charged 
Higgs-boson masses are given in GeV. 



taii/3 




ClBA 


Aew 


AsMQCD 


AsUSYQCD 


^39 


all 


5 


200 


11.241(1) 


-1.0383(3) 


-2.012(3) 


-0.00821(1) 


13.194(1) 


21.377(3) 


10 


200 


7.2568(9) 


-0.1989(5) 


-0.178(1) 


-0.00721(2) 


7.9428(5) 


14.815(2) 


20 


200 


12.546(2) 


0.1881(6) 


0.752(3) 


-0.03570(6) 


7.9968(6) 


21.447(4) 


10 


150 


12.497(1) 


-0.2574(5) 


-0.561(2) 


0.00191(4) 


8.7064(5) 


20.387(3) 


10 


400 


1.2907(2) 


-0.00530(7) 


0.0328(2) 


-0.008954(7) 


4.4386(3) 


5.7477(4) 


10 


600 


0.35740(5) 


-0.00832(2) 


0.01594(5) 


-0.006263(4) 


2.7481(1) 


3.1069(2) 



5.6 pp l^^if^: total results at 7 TeV and 14 TeV 



The total production cross section for the W final state at the LHC is shown in Fig. 11 



and Fig. 12, as well as in Table [TJ The cross section increases by an order of magnitude 



when the center-of-mass energy goes from 7 TeV to 14 TeV. The gg contribution is 
largest for small tan /3 and large Mu± while the bb dominates when tan /3 > 12 and, 



approximately, Mh± < 200 GeV. In the right panel of Fig. [TTJ one can see a little bump 
on the gg contribution around Mh± = 200GeV, attributed to the three-point Landau 
singularities discussed in Section |3| The total cross section depends strongly on the 



phases (j)t and 03 as can be seen from Fig. 12 The gg contribution is almost independent 
of 03 since the gluino does not appear at the one-loop level (the contribution through A;, 
resummation is of higher-order effect). 



The CP violating asymmetry is shown in Fig. 13 as a function of tan (3 and M^/i , and in 



Fig. 14 versus (pt and 03. The uncertainty bands obtained by varying the renormalization 



and factorization scales (we set Hr = /dp for simplicity) in the range /x_fo/2 < /i^ < 2/ii?o 



are shown only in Fig. 13 since the uncertainty depends strongly on tan (3 and in particular 
on Mh±, but not on the phases. A more detailed account of the scale uncertainty of 
our results is given in the next section. As discussed at the end of Section |4| the CP 
violating effect is dominantly generated by the gluon-gluon fusion channel. The bb channel 
contributes significantly to the symmetric cross section and thus to the denominator of 
the CP asymmetry. It is therefore easy to understand why 5cp is small for large tan P and 



small Mh±, as seen in Fig. 13 The dependence on 03 is explained by the same reasons: 
the numerator is independent of 03 while the denominator including (Tf,5 has a minimum 
at 03 = 0. The CP asymmetry is therefore maximum around 03 = 0. 

5.7 Scale dependence 

In this section we discuss the scale dependence of the total cross sections and CP asymme- 
tries. Since the calculation of the loop-induced subprocess gg — )■ W^H"^ includes only the 
leading order contribution (with improvements on the bottom-Higgs couplings and neutral 
Higgs-mixing propagators), there is no cancellation of the renormalization/factorization- 
scale dependence in this channel. We therefore concentrate on the scale dependence of 
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Figure 11: The cross section as a function o/tan/3 (left) and Mh± (right). 




Figure 12: The cross section as a function of (pt (l^ft) and 03 (right). 
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Figure 13: CP asymmetry as a function o/tan/3 (left) and Mh± (right). Within the band, 
the scale fiR = is varied in the range /ii?o/2 < Hf < '^I^fo- 
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Figure 15: The cross section (left) and CP asymmetry (right) as functions of the renor- 
malization and factorization scales (fiR = fip)- 



Table 2: Cross sections in fb for pp/bb — )■ W~H~^ and pp/gg — )■ W~H^ at different 
values of the factorization(renormalization) scale. The CP asymmetries in percentage are 
also shown. 









7 TeV 






\/s = 


14 TcV 




f^R = ^J■F 


flBA 


bb 

"NLO 




<5cp 


ClBA 


bb 

'^NLO 




Sep 


A<Fo/2 


1.1028(2) 
1.1544(1) 
1.1790(1) 


1.0434(3) 
1.0870(2) 
1.1445(2) 


1.42088(9) 
1.02168(6) 
0.7631(5) 


8.207(8) 
6.967(8) 
5.868(7) 


6.6774(8) 
7.2568(9) 
7.6648(9) 


6.633(2) 
6.873(1) 
7.224(1) 


10.4606(6) 
7.9428(5) 
6.2204(4) 


8.380(7) 
7.457(8) 
6.591(8) 



the bb — )■ W^H"^ cross section calculated at NLO, see Fig. 15 (left). We set /ir = /ip 
for simplicity. The remaining uncertainty of the NLO scale dependence is approximately 
9% (9%) when fip is varied between fiFo/2 and 2fipQ, compared to approximately 14% 
(7%) for the IBA, at 14 TeV (7 TeV) center-of-mass energy. The uncertainty is defined 
as (5 = [\a{fiFo/2) — cr{fJ'Fo)\ + Wi^ftpo) — (T{fipo)\]/a{fiFo). The IBA scale dependence 
looks quite small because we have set both renormalization and factorization scales equal, 
leading to an "accidental" cancellation. The IBA cross section increases as fip increases 
while it decreases as ^r increases. We recall that fip enters via the bottom-distribution 
functions and fiR appears in the running 6-quark mass. That accidental cancellation 
depends strongly on the value of tan/3. We have verified, by studying separately the 
renormalization and factorization scale dependence, that including NLO corrections does 
reduce significantly each scale dependence. 

Concerning the CP asymmetries, the scale dependence is shown in Fig. [15] (right). 
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We again set here = fip. If is varied between /iFo/2 and 2fiFo, the uncertainty 
is approximately 24% (34%) for 14 TeV (7 TeV) center-of-mass energy. This uncertainty 
is so large because the dominant contribution to the CP asymmetries (the subprocess 
gg — i- W^H"^) is calculated only at LO. 

In Table |2] we show the values of the cross sections for the two subprocesses as well as 
the CP asymmetries. The scale-dependence uncertainty of the gg — W^H^ process is 
indeed very large. It is mainly due to the running strong coupling as{iiR) which depends 
logarithmically on the renormalization scale. 

6 Conclusions 

In this paper we have studied the production of charged Higgs bosons in association with 
a W gauge boson at the LHC in the context of the complex MSSM. The NLO EW, 
SM-QCD and SUSY-QCD contributions to the bb annihilation are calculated together 
with the loop-induced gg fusion. Special care is dedicated to the use of the effective 
bottom-Higgs couplings and the neutral-Higgs boson propagator matrix. Moreover, the 
CP violating asymmetry, dominantly generated by the gg fusion parton subprocess, has 
been investigated. We have shown that the A;, and the Higgs-mixing resummations can 
have large effects on the production rates and CP asymmetry. 

Numerical results have been presented for the CPX scenario. It is shown that the 
production rate and the CP asymmetry depend strongly on tan /3, Mi{± and the phases 
(pt, 4>3- Large production rates prefer small tan/3, small Mh± and the phases (pt, (ps about 
±7r. Large CP asymmetries occur at small tan /3, for Mh± of about 250 GeV, and (pt ^ 
±7r/2 and 03 = 0. 

We have also studied the dependence of the results on the renormalization and factor- 
ization scales. For the bb subprocess, the NLO corrections reduce significantly the scale 
dependence while the gg fusion suffers from large scale uncertainty mainly due to the 
running as{f^R)- This makes the final results, in particular the CP asymmetry, depend 
significantly on the scales. A two-loop calculation would be needed to reduce this uncer- 
tainty to the level of a few percents. 
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A Feynman diagrams 



Wc present here the classes of Feynman diagrams that contribute to the two subprocesses 
bb W^H"^ and gg W^H"^. We use hi {i — 1,2,3) to denote the neutral Higgs 
bosons {h, H, A). 




Figure 16: One-loop SM-QCD diagrams for the partonic process bb W^H"^. 




(c) 

Figure 17: Gluon-radiation and gluon-induced QCD diagrams. 




Figure 18: One-loop SUSY QCD diagrams for the partonic process bb — >■ W^H^. 
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Figure 19: One-loop EW contributions for the partonic process bb — W^H^ . The shaded 
regions are the one-particle irreducible vertices. The diagrams are drawn with Jaxodraw 
161^. 




Figure 20: One-loop Feynman diagrams for the partonic process gg — )■ W^H^ . 
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Figure 21: Photon-radiation and photon-induced EW diagrams. 

B Counterterms and renormalization constants 

In this section, we list the Feynman rules and counterterms for vertices and propagators 
which appear in the hb — )■ W^H^ channel. They can be expressed in terms of coupling and 
field renormalization constants (RC) which relate the bare and renormalized quantities. 
The RCs are defined as in Ref. [62j for the SM-like fields and as in Ref. f27| for the Higgs 
sector. The following one-loop Feynman rules use the standard convention and notation of 



FeynArts 45 . In the vertices all momenta are considered as incoming. We introduce the 



shorthand notation Sa = sin a, = cos a, ta = tan a, sp = sin f3, cp = cos f3, tp = tan f3. 
Fermion-Fermion- S calar : 
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'HH 



'w &sw I 1 A 7 



-24^(1 + ^^e + ^ + 4^tan/3-S| 
.^(1 + + ^ + 4aan/3 - S| - ^ + 



Ssw 
Sw 



2M^ 
^kj^Z^^j^ 

4aan/3-^ 



C+ = 



Fermion-Fermion- Vector: 
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djUiW : 



k 



1 

\plsw 



[V,j{l + 5Z, 



sw 2 

5sw 1 
— — + - 
Sw 2 



Scalar-Scalar- Vector: 



S2,P2 



ieC{pi -P2y 



hH-W^ : 



hH+W- : 



HH-W^ 



HH+W~ 



AH^W^ : 



C 



C 



c 



c 
c 



cos(/3-a)r ^^^^ , 

[1 + OZe h -OZivw 



2s 



w 



Sw 



sm(/3 — a) 
2 cos(/3 — a) 



(5Z 



Hh 



SZn- 



COs(/3-a) p ^^^^1^X7 ^ ^ 

[1 + dZe h -oZww + ^OZhh + -oZh+h- 

ZSw Sw ^ z / 



sin(/3 — a) 
2 cos(/3 — a) 



sin(^ 5^ 1 

2sw ^ — 



Sw 

cos(/3 — a) 
2 sin(/3 — a 



1 1 



{SZhH + SZg-h+ 



^ [1 + oZe h -o^vK^y + ^oZhh + ^oZh+h- 

Zsw Sw z z z 



cos(/3 — a) 
2 sin(/3 — a) 



- — [l + 6Ze - + -6Zww + 7:SZaa + ^SZff+H-] 
ZSw Sw ^ z z 



4:SW 



{6Z 



AG 



5Z, 



G-H+ , 



Vector- Vector-Scalar: 




C 
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The vertices VVV^H^ with V = 7, Z do not appear at tree level. The counterterms 
are generated at one-loop level, however. 

-fW^H^ : C = ^ {SZg-h+ + sin 2(35 tan /3) , 

ZW^H^: C = -^^^^{6ZG-H++sm2/36tan(3). (43) 

2cw 

One also needs counterterms for the renormalized propagators. The complete set of 



counterterms for the scalar-scalar case can be found in Ref. 27 . We list here extra pieces 
needed in our calculation. 
Scalar- Vector: 



iMw 



{5Zg-h+ + sin 2/3(5 tan j3)p^ 



Fermion-Fermion: 



The renormalization of the fermion fields in the presence of CP violation is a bit more 
involved than the CP-conserving case. We therefore give here explicit formulae for mass 
and wave function RCs of the quark fields. The quark self-energy can be decomposed as 

= PPl^.Ap') + t^PR^^Ap') + Pl^.Ap") + Pr^.Ap")- (44) 

We note that Sg^;(p^) = T^q^rip^) = ^q,s{.P^) in the case of CP invariance. The renormalized 
self-energy is written as 



S,(p)=S,(p) + S,(p), (45) 
1 

■2' 



S,(p) =li6Z,,L + 6Z;Ai>PL + li6Z,,n + ^Zl^i^P^ 

+ ^(2^ + + 6Z;AP, + ^(2^ + 6Zl, + 6Z,n)Pn. (46) 



In general, SZg^L and SZq^R are complex[^ 6mq can always be made real by rephasing 
the field ipi (or ipR). At this step any phases can be absorbed into the two factors SZg^L 
and SZq^n which will have to be determined. It is obvious that the squared amplitude is 
invariant under a global rephasing 

il' = ipL + ipR^ e'^n^L + ^i?). (47) 

From this freedom we can, for example, make SZq^^ real while SZq^i remains complex. 
We therefore need four conditions to determine the three renormalisation constants. The 

^If we impose CP invariance then SZq^L and SZq^n can be taken real. 
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OS conditions are 

Retq{p)u{mq) = 0, 



j6 — m, 



Ref,(p) 



UiTTlr, 



iu{mg), 



Reti(mq)fg(p) = 0, u{mg 



p-rriq 



iu[m, 



(48) 
(49) 



where Tq{p) = — ra^ + Sg(p)] and = m^. Re sets the imaginary part of the loop 
integrals to zero since they are not involved in the renormalisation. The Hermiticity of 



the Lagrangian imposes 63 



^g{p) = -7(lf5(p)7o- 



(50) 



It is obvious that Eq. (49) can be derived from Eq. (48) and Eq. (50). From these 
conditions we get the following results 

Srug = ^Re|mg[Sg,i 
SZg^L '■ 



[m 



-'q,T\"<'ql f ) 



1 



-Re^S,,(mJ)- — 



(Sg,z(mJ) - Sg,^(mJ)) 



+ m, 



d 



K(S,,l(p') + S,,«(p2)) + E,Xp') + S,,,(p2)]^,^^J, 



Re I Y:q^R{ml) 



2m, 



■(S,,^(mJ) -Sg/mJ)) 



+ m, 



"dp^ 



(51) 



where we have used the freedom Eq. (47) to take lm{5Zg fi) = — lm(5Zg^L). These results 



agree with the ones in 21 
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